A Novel Test of the Modified Newtonian Dynamics with Gas Rich Galaxies 
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The current cosmological paradigm, ACDM, requires that the mass-energy of the universe be 
dominated by invisible components: dark matter and dark energy. An alternative to these dark 
components is that the law of gravity be modified on the relevant scales. A test of these ideas is 
provided by the Baryonic TuUy-Fisher Relation (BTFR), an empirical relation between the observed 
mass of a galaxy and its rotation velocity. Here I report a test using gas rich galaxies for which 
both axes of the BTFR can be measured independently of the theories being tested and without the 
systematic uncertainty in stellar mass that affects the same test with star dominated spirals. The 
data fall precisely where predicted a prion by the modified Newtonian dynamics (MOND). The 
scatter in the BTFR is attributable entirely to observational uncertainty. This is consistent with 
the action of a single effective force law but poses a serious fine-tuning problem for ACDM. 



The mass discrepancy problem in extragalactic sys- 
tems is well established. When known dynamical laws are 
applied to these systems, the observed mass in stars and 
gas falls well short of explaining the observed motions. A 
classic example is that the rotation curves of disk galax- 
ies tend to become roughly flat {Vf ^ constant) when 
they should be falling in a Keplerian {V cx r~^/^) fash- 
ion. The common interpretation for this phenomenon is 
dark matter. However, a logical alternative is that the 
dynamical laws that lead to the inference of dark matter 
need to be revised on the scales appropriate to galaxies. 

One striking fact about extragalactic systems is that 
they are many orders of magnitude larger than the solar 
system in which conventional dynamics is extraordinarily 
well tested. One idea is thus to modify gravity on some 
suitably large length scale such that the apparent need 
for dark matter would be manifest in galaxies but not 
in the solar system. Such size-dependent ideas fail and 
can generically be excluded However, there are other 
ways in which galaxies differ from the solar system. For 
example, the centripetal acceleration required to keep a 
star in orbit in a galaxy is very much lower than that ex- 
perienced by the planets orbiting the sun: ^ 10^^° ms^^ 
vs. 6 X 10"'^ ms~^ for the Earth. 

MOND posits a new constant with dimensions of ac- 
celeration, ao, which defines the boundary between con- 
ventional dynamics and a new domain of dynamics. The 
conventional dynamics hold in the limit of high accelera- 
tion, a ^ flo, and the modified regime occurs in the limit 
of low accelerations, a ^ an. The value of ao must be 
determined observationally [3|, but once specified is con- 
stant. In the modified regime, rotation curves become 
asymptotically flat far from a central mass p. This fol- 
lows from the scale invariance symmetry of the equations 
of motion under transformations {t,r) — > {\t, \r) 4]. An 
absolute relation between the asymptotically flat rotation 
velocity Vf and the total mass 



aoGMb = vf 
follows uniquely on dimensional grounds. 



(1) 



Rotationally supported galaxies follow an empirical re- 
lation between mass and rotation velocity known as the 
Baryonic Tully-Fisher Relation (BTFR) fs']. This empir- 
ical relation can, in principle, provide a quantitative test 
of the prediction of MOND. In order to do so, we require 
independent, accurate measurements of both Mt, and Vf. 
While the latter is readily obtained from resolved rota- 
tion curves, mass determinations are more problematic. 

The baryonic mass is the sum of both stars and gas: 
Mb — Mi, + Mg. A great deal is known about stars, but 
stellar mass estimates for galaxies are subject to a sys- 
tematic uncertainty of ^ 0.15 dex because of uncertainty 
in the stellar mass function and some details of stellar 
evolution This level of systematic uncertainty pre- 
cludes an unambiguous test of ([1]) with star dominated 
spirals [I]. 

A clean test of the BTFR predicted by MOND follows 
if we can identify a class of galaxies where stars do not 
dominate the baryonic mass budget. Atomic gas typi- 
cally dominates the mass of non-stellar material in disk 
galaxies. Its mass follows directly from the distance to 
each galaxy, the measured 21 cm flux, and the physics 
of the spin flip transition of hydrogen. It does not suffer 
from the systematic uncertainty of stellar mass. 

Late type, low surface brightness disk galaxies fre- 
quently have gas masses in excess of their stellar masses 
[Sj. When Mg > M^,, the systematic uncertainty in stel- 
lar mass is reduced to a minor contributor to the er- 
ror budget (Fig. [ij. Thanks to recent work it is 
now possible to assemble a large sample (47) of galaxies 
with Vf measured from resolved rotation curves that sat- 
isfy the gas domination criterion Mg > Mi,. This prop- 
erty enables a novel test of MOND with no free parame- 
ters. Both Mb and Vf are directly measured and, for the 
first time, are not dominated by systematic uncertainties. 
Moreover, these galaxies are unambiguously in the deep 
modified regime where ([J) holds, with Vj /vmax ~ ao/10. 
Here the distinction between MOND and ACDM is most 
pronounced. 

Fig. [2] shows the gas rich galaxy data together with 
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FIG. 1. (a) The mass of stars and gas in rotating galaxies. Triangles represent star dominated spirals [3| with > Mg. The 
data for gas rich galaxies with Af* < Mg come from several independent sources denoted by circles [9|, squares [l3|, and stars 
111 ], (b) The fraction of the error budget asys/utot contributed by the systematic uncertainty in stellar mass {cTsys ~ 0.15 dex 
g|) as a function of the total baryonic mass Mb — M* + A'Ig. Our knowledge of the masses of star dominated galaxies is limited 
by this systematic uncertainty, but it has little effect on gas dominated galaxies. 



the predictions of MOND and ACDM. The data fah pre- 
cisely where MOND predicts. This happens with no 
fitting whatsoever — there are zero free parameters in 
Fig. [2] Computing with the slope fixed to 4 and the 
normahzation fixed at the previously determined value 
of ao 01 gives = 44.3 for 46 degrees of freedom for 
a reduced xl — 0.96. If we treat Oq as a fit parame- 
ter and minimize with the slope fixed to 4, we find 
ao = 1.24 ± 0.14 X 10~^° ms~^. This is indistinguish- 
able from the previous value, and actually increases 
because we have added an unneeded degree of freedom: 
xl = 0.99. If we further treat the slope as an additional 
free parameter, we find 3.8 ± 0.2. This does not differ 
significantly from the MOND prediction of 4, nor does it 
improve the fit: xt — 0.98. The data therefore provide 
no reason to suspect a BTFR that differs in any way from 
that predicted by MOND. 

The specific BTFR that the data follow is unique to 
MOND. Indeed, to the best of my knowledge, MOND is 
the only theory to make a strong a 'priori prediction for 
the BTFR. The dark matter paradigm makes no compa- 
rably iron-clad prediction. 

The expectation in ACDM is that total mass (both 
dark and baryonic) scales with rotation velocity as Ma = 
(A/2)-i/2(G'i^Q)-i ]/|. These quantities are defined at 
a radius where the enclosed density exceeds the cosmic 
critical density by a factor A. The virial radius occurs 
at A « 100 [12]. This notional radius is well beyond the 
reach of observations. To plot the ACDM line in Fig. [21 
we assume V/ = Vvir and Mf, = h^vir where /(, = 0.17 
is the cosmic baryon fraction [l3|. This nominal expec- 



tation has the wrong slope and the wrong normalization. 

In order to reconcile ACDM with the data, we must 
invoke additional parameters. The simplest assumption 
is that only a fraction of the baryons in a halo are 
detected: Af^ — fdftMyir. Once we have granted our- 
selves this freedom, a galaxy could, in principle, have 
any fd<i and reside anywhere below the ACDM line in 
Fig. [21 From this perspective, it is puzzling that galaxies 
reside only along the line predicted by MOND. 

Reproducing the observed BTFR in ACDM requires 
a remarkable degree of fine-tuning. The detected 
baryon fraction must follow the formula log = 
log (V// 100 kms^^) — 1.2. Astrophysical feedback is of- 
ten invoked to cause fd < I, but provides no satisfactory 
explanation for why this particular tuning of fd arises. 

A further test is provided by the scatter in the observed 
relation. In ACDM, any scatter in fd translates directly 
into the BTFR: there should be at least some intrinsic 
scatter. In MOND, the BTFR is a consequence of the 
force law, and should have no intrinsic scatter. 

Fig. [3] shows a histogram of the ratio a = V^/{GMh) 
formed from the data. This represents the scatter around 
the MOND line in Fig. [21 If the data are randomly 
distributed, they should approximate a gaussian whose 
width is dictated by the size of the error bars. Such a 
gaussian is shown in Fig. [3l It is not fit to the data; it 
is simply centered on the previously determined value of 
ao [3] with a width corresponding to the uncertainty in 
the data. Only random errors are considered here; the 
residual systematic uncertainty in the stellar mass corre- 
sponds to a small shift in the total mass and should not 
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FIG. 2. The BTFR for gas dominated galaxies. The sum of 
detected baryonic mass, stars and gas, is plotted against the 
fiat rotation velocity V/ (symbols as per Fig. [ij. Both mass 
and velocity are measured independently of either MOND or 
ACDM. The data are well removed from the expectation of 
the standard cosmology (upper line) , but follow the prediction 
of MOND (lower line) with no fitting whatsoever. 



introduce additional scatter. 

Observational uncertainty suffices to explain the scat- 
ter in the data. The data are consistent with a BTFR of 
zero intrinsic width. This is natural if the BTFR is im- 
posed by the force law, as in MOND. It is not expected 
in ACDM where there should be many sources of scatter. 

From the perspective of cosmology, it is disturbing that 
MOND works at all. If ACDM is the correct paradigm, 
this should not happen 



14| . Yet when pressed into a 



new regime where the predictions of the two theories are 
distinct, MOND outperforms ACDM. 

This is not the first time that strong predictions of 
MOND have been realized. For example, MOND pre- 
dicted in advance that galaxies of both high and low sur- 
face brightness would fall on the same BTFR 15, [l6| . 
contrary to the natural expectation of purely Newtonian 
gravity [ll[lil. It is well established that MOND pro- 
vides good fits to the detailed shapes of rotation curves 
with only the stellar mass-to-light ratio as a free param- 
eter [l9| . The required mass-to- light ratios are in good 
agreement with stellar population synthesis models [20|. 
A simple model motivated by MOND provided the only 
successful a priori prediction of the first-to-second peak 
amplitude ratio of the acoustic peaks of the cosmic back- 
ground radiation: Ai.,2 = 2.4 predicted f^l*! vs. 2.34±0.09 
measured [22j. It is rare for a non-canonical theory to 
have so many predictive successes. 

MOND also has its share of problems. The same ansatz 
that correctly predicted the second acoustic peak ampli- 



tude also predicts a lower third peak than is observed 
This does not falsify MOND, but it does imply that 
a generally covariant parent theory should provide an ef- 
fective forcing term [2j| . 

The most serious observational problem facing MOND 
is the dynamics of rich clusters of galaxies. These appear 
to weigh more than can be accounted for with the ob- 
served baryons even with the modified dynamics 25, 26| . 
This residual mass discrepancy is roughly a factor of two 
in mass. On the one hand, this is very disturbing — a 
theory that seeks to eliminate the need for cosmic dark 
matter itself suffers a missing mass problem. On the 
other hand, this is less severe than the missing baryon 
problem in ACDM, where dwarf galaxies are missing 
99% of the baryons that should be associated with their 
dark matter halos So both theories suffer a missing 
baryon problem, albeit of different amplitudes in systems 
of vastly different scale. 

While some of the mass in clusters appears to be dark, 
even in MOND, there is nothing that requires this unseen 
mass to be in some new form of non-baryonic particle. 
Indeed, big bang nucleosynthesis implies the existence of 
considerably more baryons than have so far been detected 
[iit . If only a fraction of these missing baryons reside 
in clusters it would suffice to resolve the residual mass 
discrepancy suffered by MOND. 

Perhaps the most prominent example of a cluster with 
a serious residual discrepancy in MOND is the bullet clus- 
ter [2^ . In this system, the gravitational lensing of back- 
ground galaxies indicates that the mass is offset from the 
X-ray plasma. This is the same residual mass discrep- 
ancy that is seen in all rich clusters. While the bullet 
cluster is frequently cited as evidence against MOND, 
it is also problematic for ACDM. The sub-clusters that 
compose the bullet cluster collided at a remarkably high 
velocity (~ 4700 kms~^). This is exceedingly unlikely in 
ACDM, occurring with a probability of only a few parts 
in a billion [s^. In contrast, such high collision veloci- 
ties are natural to MOND 3l[. Taken at face value, the 



bullet cluster would seem to simultaneously support and 
falsify both theories with equal vigor. 

Given the nature of astronomical data, some excep- 
tions to any theory are to be expected. What is surpris- 
ing in the case of MOND is that it continues to enjoy 
predictive successes at all. These motivate the search 
for a more complete gravitational theory that contains 
MOND in the appropriate limit 32l435|. 

It is possible that non-baryonic cold dark matter does 
not exist. If it does, and ACDM is the correct solution, 
the challenge is to understand the empirical systematics 
encapsulated in the simple MOND formula. These are 



not native to the current cosmological paradigm 1J| but 
must be explained by any successful theory. 

Another possibility is that dark matter particles have 
properties that impose MOND-like phenomenology 
13811 . In this case, it is desirable to have dark matter that 
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FIG. 3. Histogram of the measured values a — Vf /{GMi,). 
The data are consistent with a normal distribution (smooth 
curve) that is centered on the previous determination of ao 
[3| with a width specified by the mean uncertainty a = 0.24 
dex. This consistency imphes a universal acceleration scale 
with negligible intrinsic scatter. This is expected in MOND, 
but poses a fine-tuning problem for ACDM. 



behaves like standard cold dark matter on large scales, 
but which interacts with normal matter so as to impose 
the MOND phenomenology in galaxies. This suggests 
some strong new form of interaction between dark matter 
and baryons. 

If MOND is essentially correct in that it is pointing 
towards an extension of gravitational theory, then the 
experiments seeking to detect dark matter will find null 
results. This would also be the case if dark matter has a 
different nature than currently presumed. If dark matter 
is detected, then the MOND formula is still useful as 
a phenomenological constraint on the effective force law 
in spiral galaxies. In any case, the predictive power of 
the simple formula proposed by Milgrom [2] is telling us 
something profound about Nature. 
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grom and Bob Sanders for sharing their perspectives, and 
the referees for their constructive comments. The work 
of SSM is supported in part by NSF grant AST0908370. 
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